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Abstract 

(N | 

We discuss the factorization breaking observed in diffractive dijet photoproduction by 
the HI and ZEUS collaborations at HERA. By considering the effects of rapidity gap 



survival, hadronisation, migration and NLO contributions, we find that the observed data 
are compatible with theoretical expectations. 
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1 Introduction 



As is well known, in QCD the cross section for a 'hard' inclusive process factorizes into universal 
parton densities and calculable hard subprocess cross sections. However, for diffractive processes 
the factorization into universal diffractive parton densities and the known subprocess cross 
sections may be broken, since the rapidity gaps associated with the diffractive process can be 
populated by secondary particles from 'soft' rescattering. 

Here, we discuss factorization breaking of the diffractive photoproduction of dijets at HERA0, 
where each jet has large transverse momentum, Ep. That is, events with a large rapidity gap 
between the proton and the hadronic dijet system produced by a photon with virtuality Q 2 ~ 
(see for example [TJ Hj). The domains xp < 0.03 and xp < 0.025 were selected by the HI 
and ZEUS collaborations, respectively. The comparison of these data with theory is not well 



x The HERA data for the difFractive production of dijets in deep inelastic scattering are consistent with NLO 
predictions [1, 2, 3]- That is, within the present uncertainties of the data and theory, no evidence of factorization 
breaking is observed. 



understood. The situation may be summarised as follows. In general, there is a tendency for 
the observed cross sections to be smaller than those predicted [51 El E] by NLO QCD. Indeed, 
for the HI choice of jet cuts, -Erjeti > 5 GeV (-&rjct2 > 4 GeV), the ratio of data/theory is about 
0.5 — 0.6 independent of the observed x 7 p], indicating an overall suppression relative to the 
QCD prediction. On the other hand, with higher jet cuts, -Erjeti > 7.5 GeV (-&rjet2 > 6.5 GeV), 
the data of the ZEUS collaboration [I] give a data/theory ratio of 0.9 if the diffractive PDFs of 
HI fit B [5] are used, compatible with little or no overall suppression^]. Moreover, in the latest 
HI analysis [10] the events have been selected using a similar set of cuts to those adopted by 
ZEUS, with, in fact, identical choices of the Et and xp cuts. In this case, the (preliminary) HI 
results give, using the HI 'jets' diffractive PDFs [2^, a data/theory ratio consistent with an 
overall suppression of about 0.8 [TO], which, within the 20 — 30% experimental uncertainties, is 
not in contradiction with the findings of ZEUS. 

There was an attempt to describe the factorization breaking by soft spectators from the 
photon interacting inelastically with the proton target and producing secondaries which popu- 
late the rapidity gap. For the hadron-like component of the photon wave function this would 
have produced a suppression by a factor of about 3, corresponding to a gap survival factor 
S 2 = 0.34 [TTj . This idea was widely discussed and used, for example, in the studies of Klasen 
and Kramer [HIE]. Nevertheless, the absorption of the hadron-like component of the photon 
cannot explain the suppression of the dijet yield observed in the largest x 7 bin, close to x 1 — 1. 
When the dijet system carries away almost all of the incoming photon momentum, the 'direct 
7 — > dijet' subprocess dominates and we anticipate that S 2 ~ 1, since the absorptive cross 
section of the point-like direct photon is very small. 



2 Components of the photon 

Before continuing the discussion, it is useful to elaborate the different components of the pho- 
ton, since they are partly a matter of convention. First, we have the 'direct' and 'resolved' 
contributions of photon interactions. However, recall that the PDFs which describe the par- 
ton distributions of the 'resolved' photon contain not only the 'hadron-like' component, but 
also a so-called 'point-like' contribution which originates from the inhomogeneous term in the 
DGLAP evolution equation for the photon [12]. From a Feynman-graph viewpoint, this 'point- 
like' contribution looks just like the diagrams for the 'direct' photon interaction. However, the 
generally accepted convention is to include the part of this contribution which corresponds to 
virtualities, \q 2 \, lower than the factorisation scale fi 2 F , in the PDFs of the 'resolved' photon, 

2 In a recent ZEUS paper [9] the data on diffractive dijet production in DIS have been included in a DGLAP 
analysis to better constrain the diffractive gluon densities. With these diffractive PDFs the dijet data are well 
described without factorization breaking, but see the comments in the concluding Section. 

3 If HI fit B diffractive PDFs are used then the ratio for the preliminary HI data implies an overall suppression 
of about 0.6 [TO]. 
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while the part corresponding to larger virtualities, \q 2 \ > fi 2 F , is included in the NLO matrix 
element of the 'direct' photon interaction. 

Unlike the 'hadron-like' component of the resolved photon, the probability of an additional 
soft interaction for the 'point-like' component is rather small. It is driven by the size of the qq 
dipole produced by the initial photon, 7 — > qq. Thus the absorptive effect for the 'point-like' 
component of the photon PDFs is not so strong; we expect the gap survival factor S 2 to be 
rather close to 1. The fact that the 'point-like' component of the resolved photon has a different 
absorptive cross section to that for the 'hadron-like' part, and therefore that the 'point-like' 
component should have a different gap survival factor S 2 , was discussed first in [5]. 

Moreover, a detailed study of the photon PDFs, that have been used to make predictions, 
show that practically in the whole of the observed x 1 region selected by the HI and ZEUS cuts 
(that is, cuts on pTjet and r/j e t, which lead to a; 7 > 0.1 — 0.2), the contribution which originates 
from the point-like photon dominates (see Fig. 4 of GRV [H]). As mentioned above, this 
component of the photon wave function has no soft spectators and the only absorption comes 
from the rescattering of the relatively small-size qq dipole that is directly produced by the 
7 — > qq transition via the inhomogeneous term in the DGLAP evolution of the photon PDFs. 
To account for the absorptive effect, this inhomogeneous contribution should be multiplied by 
the well-known eikonal factor exp(— <j a b s /-R 2 ), which is of analogous form to that introduced 



in the saturation model of [2]. Here, the absorption cross section, cr a b s , depends on the size 
of the pair, which, in turn, is driven by the running scale fi in the DGLAP evolution, that 
is cr a bs oc 1/yU 2 . In other words, in diffractive processes the point-like photon contribution 



should be suppressed by S ~ exp(— a/fx ) inside the DGLAP evolution, where the value of 
the parameter a, namely a ~ 0.6 GeV 2 , has been obtained from the known u q q- p of the dipole 
mo 



Clearly this suppression is only non-negligible near the starting scale of the evolution, 



So to summarize: 

• for the hadron-like component we expect S 2 ~ 0.34, 

• for the direct photon contribution we have S 2 = 1, 

• for the point-like component of the resolved photon, generated by the inhomogeneous 
DGLAP term, S 2 should be accounted in the DGLAP evolution. That is, it will be of 
'power-correction' form, S 2 ~ 1 — a/fi 2 + which will only be non- negligible at the 
beginning of the evolution. 

In the existing photon PDF analyses this power correction is neglected. In the more recent 
'AFG' photon PDF analysis [18], where the optimum starting scale was /ig = 0.7 GeV 2 , the 
corresponding effect is practically negligible (< 20%). On the other hand, we expect that 
the use of the 'GRV [13J photon PDFs, which have a starting scale of = 0.3 GeV 2 , will 
overestimate the point-like part of the resolved photon contribution (by up to about 25%). 



4 This value, a ~ m' ], is consistent with the dipole parametrizations of Refs. [T51 [TBI H7) . 
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3 Interpretation of HERA data 



We return to the discussion of the measurements of the diffractive photoproduction of dijets 
at HERA that we outlined in the Introduction. After all the cuts have been made, the point- 
like component of the photon dominates (except in the largest x 7 bin, where 70% comes from 
the direct component^) . We therefore expect only weak absorption throughout the observed x 7 
interval. Indeed, this expectation of only a small amount of factorization breaking for the higher 
E T dijet data (with -Erjeti > 7.5 GeV and E T - ]et2 > 6.5 GeV) data is, within the uncertainties, 
compatible with the observations of ZEUS [I] and with an overall suppression of about 0.8 



data/theory to be larger in the largest a; 7 bin, though the variation is within the uncertainties 
of the ratio. 

Most probably, in the largest x 7 bin, we observe the interplay of a few different physical 
effects. First, the GRV photon PDFs [32] used in the theoretical calculations, overestimate the 
effect of the point-like contribution of the resolved photon due to a too small starting scale Hq. 
Second, as described above, there should be some absorption in the inhomogeneous DGLAP 
term at the beginning of the evolution. Even using the AFG photon PDFs [18J, we expect 
a bit of absorption. From these two effects, we expect some weak factorization breaking in 
the region where the resolved photon contribution dominates. However, in the largest x 7 bin, 
we might anticipate no suppression for the direct photon component. On the other hand, the 
hadronization corrections estimated from the LO RAPGAP Monte Carlo [19] . and applied to 
the NLO theoretical calculations of Refs. [5j EJ [7], are rather large. In comparison with the 
value of x 7 calculated at the parton level, after hadronization the fraction x 7 of the photon's 
momentum carried by the dijet system becomes smaller. Events from the largest x 7 bin migrate 
to lower Xj. Therefore, we expect a negative correction due hadronisation in the largest x 7 bin 
and a positive correction for x 1 ~ 0.7. Thus, the final effect brings the event rate in the largest 
x 7 bin into line with the other x 7 bins, providing, surprisingly, an almost uniform suppression 
independent of x 7 for x 7 > 0.2. 

Recall that for large photon virtualities, that is, for diffractive high-.&r dijet production 
in deep inelastic scattering (DIS), the data [TJ [2J [3] are consistent, within uncertainties, with 
the NLO predictions without any factorisation breaking. Here we deal with a very small- 
size incoming qq dipole and so the absorptive effects are negligibl^]. For this DIS case, it 

5 This can be estimated in various ways - for example, from the results of Monte Carlos or by inspection of 
the results in [5]. 

6 Recall that we noted that the ZEUS [4] and preliminary HI data [10] were obtained using HI fit B [8] and 
HI 'dijet' [5] diffractive PDFs respectively. 

7 Note that there is no suppression and, correspondingly, no factorisation breaking, in the case of diffractive 
charm photoproduction where: (a) the initial heavy-quark small-size dipole, cc, has a small absorptive cross 
section, and (b) the value of x 7 is not fixed - so hadronisation does not affect the results. Next, recall that the 
ZEUS data [4], obtained with the larger Et threshold, are consistent with no, or little, suppression in the whole 
observed range of x 1 . 




Moreover, there is a (weak) tendency of the ratio of 
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would be interesting to study the role of the hadronisation and migration as a function of x 1 . 
Unfortunately, in [TJ [2] only cross sections integrated over x 7 are presented, so the effects of 
migration are not revealed. The diffractive DIS x 1 dijet distribution was shown in [3J, see also 
[2~T] . At large x 7 the data agree with the RAPGAP LO Monte Carlo. However, the NLO 
calculation leads to a larger cross section than that at LO level. Indeed, in Fig. 6(d) of [3J the 
NLO prediction in the largest x 7 bin exceeds the data by about 20% or more. This may be 
interpreted as an indication in favour of migration (similar to that in photoproduction) which 
decreases the observed value of x 1 . 

We also note, from detailed studies made by Klasen and Kramer [22J, that at the NLO 
level, a rather large fraction of the events correspond to three-jet production. For example, the 
fraction is up to 30% if the factorisation scale is chosen to be /zp = Et and the cut on transverse 
momentum of the third jet is taken to be p^T > 1 GeV. Of course, in three-jet events, the two 
highest Et jets cannot carry x 7 = 1. On the other hand, the value, x^, of momentum fraction 
carried by the third (lowest Et) jet, cannot be too large. So these three-jet events go mainly 
to the bins with x 7 = 0.6 — 0.9. Moreover, for the large factorisation scale /zp = Et, these 
three-jet events are in an azimuthal 'Mercedes'-like configuration. This is unlike the additional 
jets generated as parton showers by RAPGAP (and other) Monte Carlos, which are collinear 
to a high Et or beam jet. Thus, the hadronisation corrections calculated based on the LO 
Monte Carlo should be different from the effects associated with NLO contribution. In the 
'Mercedes'-like NLO configuration it is natural to expect stronger migration from the highest 
x 7 bin to smaller values of x 7 . 

Note that using lower Et cuts on the jets, we expect both the hadronisation corrections 
and the absorptive effect to increase. As a rule, the non-perturbative hadronisation correction 
reduces the value of Et by some constant amount AEt- Therefore the effect of hadronisation 
is expected to be stronger for lower Et- Next, a larger part of the DGLAP evolution is affected 
by the absorption of the qq dipole for lower Et- This is because the dipole size, ~ l//z with 
/z ~ a, is not small in a larger part of the interval in In /z 2 space going from /zo up to Et, see 
eq. ([I]) and footnote 7 below. Indeed, for the lower choice of Et cuts, the data/theory ratio 
is observed to be smaller (~ 0.5 — 0.6) [TJ. However, only a minor part of this suppression is 
predicted to come from soft rescattering, S 2 . 

For a simple estimate we may use S 2 = exp(— a//z 2 ) with a = 0.6 GeV 2 to calculate the 
expected suppression of the point-like component of the resolved photon. We use the 'AFG' 
photon PDFs [18]. Then for the quark- and gluon-partons we find_| 

S 2 = 0.84(0.81), S 2 = 0.74(0.71) (1) 

for data samples with ETjeti > 7.5(5) GeV. The quark contribution dominates for large x 7 , 
and even for x 7 = 0.2 the gluon contribution is only 20%. Thus, using the value of S 2 , we 
expect the absorptive suppression to be S 2 ~ 0.8, except in the largest x 7 bin, where the direct 

8 Using the 'GRV photon PDFs [13] the corresponding predictions are S 2 q = 0.75(0.71), Sf = 0.58(0.53). 
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contribution (with S 2 = 1) dominates; the resolved contribution from the quark gives only 30% 
so the final prediction is S 2 ~ 0.94 in this bin. Then we have to consider hadronisation and 
the effects of migration. Of course it would be better to implement the NLO '2 — > 3' matrix 
element into a Monte Carlo generator (say, RAPGAP) and to trace the result of hadronisation 
of the three-jet system explicitly^. At the moment we have no such possibility. However, it is 
reasonable to assume that due to the non-negligible three-jet contribution, the hadronisation 
correction should be about 10% stronger than that given by the present RAPGAP, decreasing 
(by about 10%, as indicated by the hadronization corrections in the RAPGAP LO Monte Carlo) 
the NLO prediction for largest x 1 bin and enlarging (by about 5%) the predictions for lower 
x T This, together with the S 2 suppression of the point-like component of the resolved photon, 
can satisfactorily describe the present experimental data. 



4 Conclusions 

If we were to neglect all final state interactions, and have no hadronisation corrections or migra- 
tion effects, we would have exact factorisation in the DIS diffractive production of dijets, but 
some factorisation breaking in diffractive photo-produced dijets where the hadron-like compo- 
nent of the resolved photon suffers a small gap survival probability, for example, S 2 ~ 0.34 [llj. 
However in the kinematic region of the HERA dijet data, x 7 > 0.1, the hadron-like contribution 
is very small. Here the point-like part of the resolved photon dominates, for which S 2 is close 
to one, see (pQ). 

We have argued that the factorisation breaking observed in the diffractive photoproduction 
of dijets is most probably caused by hadronisation corrections and migration effects in the final 
state. These effects should be similar in both DIS- and photo-produced diffractive dijets. For 
this reason it is dangerous to include DIS diffractive dijet data in diffractive PDF (DPDF) 
analyses. The analyses assume exact factorisation and therefore tune the DPDFs to partly 
compensate for the effect of the interactions in the final state. Note that after including DIS 
diffractive dijet data in the DPDF analysis, ZEUS [9] do not observe factorisation violation 
in diffractive photoproduced dijets. This observation may be considered as support for our 
argument. That is, the factorisation breaking, observed when diffractive dijet photoproduction 
data are described using DPDFs obtained by analysing pure inclusive diffractive DIS data (such 
as HI fit B [8] or MRW [24]), is mainly of hadronisation/migration origin. 

In summary, the hadron-like component of the resolved photon, which is suppressed by a 
factor S 2 ~ 0.34 [TT] , only starts to be important for small x 7 . Indeed, to feel the hadron- 
like component one needs to observe dijets far in rapidity from the photon, corresponding to 
x 7 < 0.1. This region was difficult to access at HERA 10 !. The point-like component of the 

9 In the ideal case, the absorption of the incoming qq dipole in the point-like part of the resolved photon 
component, that is the factor S 2 = exp(— a//i 2 ), should be included in this Monte Carlo as well. 

10 This statement is not strictly true. For example, in Ref. [23], the HI collaboration presented results down 
to x 1 = 0.05 for dijets with jets of Et jet > 4 GeV and |77j e t| < 2.5. 
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resolved photon, which is calculable perturbatively, is the dominant one for x 7 > 0.1, and has 
a small suppression. For this component, the spectator partons have relatively large transverse 
momenta and can be seen as a third jet. Finally, after including the direct component and 
taking into account the effects of hadronisation and migration, we find that our expectations 
are consistent with the observed data for diffractively photoproduced dijets. 
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